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In this work we explore the possible evolutionary track of the neutron star 
in the newly discovered Be/X-ray binary SXP 1062, which is believed to be the 
first X-ray pulsar associated with a supernova remnant. Although no cyclotron 
feature has been detected to indicate the strength of the neutron star's magnetic 
field, we show that it may be > 10 14 G. If so SXP 1062 may belong to the 
accreting magnetars in binary systems. We attempt to reconcile the short age 
and long spin period of the pulsar taking account of different initial parameters 
and spin-down mechanisms of the neutron star. Our calculated results show that, 
to spin down to a period ~ 1000 s within 10 — 40 kyr requires efficient propeller 
mechanisms. In particular, the model for angular momentum loss under energy 
conservation seems to be ruled out. 
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1. Introduction 

K 

As a subgroup of high- mass X-ray binaries (HMXBs), Be/X-ray binaries (BeXBs) consist 
of a neutron star (NS) and a Be companion star which show emission lines and infrared (IR) 
excess in its spectrum. The origin of the emission lines and the IR excess is attributed to the 
circumstellar discs, which is fed from the material expelled from the rapidly rotating Be stars. 
The X-ray emissions is believed to originate from accretion of matter in the circumstellar 



discs by the NSs (see lReigll201ll for a review). 



BeXBs are subdivided into persistent and transient sources according to their different 
X-ray properties. Transient systems are characterized by outbursting activities, in which 
the X-ray flux increases by ~ 1 — 4 orders of magnitude compared with in quiescent state, 
and the outbursts typically last about 0.2 — 0.3 orbital period. These systems often have 
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moderately eccentric (e > 0.3) and relatively narrow orbits (P OT b ^ 100 days). Persistent 
sources are relatively quiet systems with low X-ray luminosities (Lx ~ 10 34 — 10 35 ergs _1 ), 
the variability of which is less than an order of magnitude. They usually have slowly rot ating 
NSs (P > 200s), low eccentric (e < 0.2) and relatively wide orbits (P or b ^ 200 days) (IReig 



201lh . 



Recently iHenault-Brunet et al.l (120121 ) reported a new BeXB SXP 1062 in the Wing of 
the Small Magellanic Cloud (SMC). This source was first discovered as a transient BeXB 
during the XMM-Newton and Chandra observat ions in 2010, and was not active dur i ng th e 
ROSAT and ASCA observations of the SMC JHaberl et al.l I2OO0I : lYokogawa et~aD l2003h . 
However, SXP 1062 seems to share some characteristics with persistent BeXBs: a relatively 
low intrinsic X-ray luminosity Lx — 6.3(^o'g) x 10 35 ergs _1 (corresponding an accretion 
rate of M = Lx/i]c 2 ~ 6 x 10 15 gs _1 with energy conversion efficiency i] = 0.1, c is the 
speed of light), a slowly rotating NS with period of P ~ 1062 s, a relatively flat light curve 
with sporadic fluctuation less than an order of mag nitude, and a probably wide orbital 
period P Q rb ~ 300 days derived from the ICorbetl (119841 ) diagram. What makes this discovery 
most noticeable is that SXP 1062 is located in the center of a shell-li ke nebula, which is 
considered to be a supernov a remnant (SNR), aging only ~ 10 — 40 kyr ( jHaberl et al.ll2012 
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Henault-Brunet et al.ll2012l ). Thus SXP 1062 provides the first example of an X-ray pulsar 
associated with a SNR, and it challenges the traditional spin-down model of NSs because of 
its extraordinarily long spin period combined with a relatively young age. 



Haberl et al.l ( 120121 ) measured the spin period change in SXP 1062 over a 18 day duration 
of observation. Their timing analysis shows that the NS in SXP 1062 has a very large average 
spin-down rate with the spin frequency derivative i> ~ 2.6 x 10~ 12 Hzs _1 (or period derivative 
P ~ 100 syr" 1 ). If the NS has a normal magnetic field {B ~ 10 12 - 10 13 G), it's hard to 
spin-down to a period ~ 1000 s within a few 10 4 years. Assuming in the extreme case that 
the NS has spun d own with a c o nstan t spin frequency derivative —2.6 x 10~ 12 Hzs _1 over 
it's whole lifetime, lHaberl et al.l ( 120121 ) derived the lower limit of the initial spin period of 
the NS is 0.5 s. Since the duration of the observation lasts only 18 days (probable less than 
one tenth of the orbit period), the extraordinarily large spin-down rate is very unlikely to 
sustain in the whole lifetime of SXP 1062, thus the present value of spin-down rate may be 
just a short-term one. 



Popov fc Turollal (120121 ) (hereafter PT12) suggest another possibility to reconcile the 



long spin period and short age of SXP 1062. Assuming that the NS is spinning at the 
equilibrium peri od, PT12 estimated t he current magnetic field to be B < 10 13 G according 



to the model of IShakura et al.l (120121 ). Their calculations show that if the NS in SXP 1062 
was born as a magnetar (B > 10 14 G) it can be spun down to ~ 1000 s within a few 10 4 yr. 
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In this work we consider the proposed mechanisms that can account for the observed 
rapid spin-down in SXP 1062, if it is alternatively a magnetar with current field strength 
> 10 14 G. In Section 3 we point out that this possibility remains according to current 
observations. Based on this, we investigate its spin-down evolution, taking account of various 
kinds of braking torques during the propeller stage, which is introduced in Section 2, to 
examine how the observations of SXP 1062 can constrain the possible spin-down mechanisms 
in NSs. We present our calculated results in Section 4, and discuss their possible implications 
in Section 5. We conclude that the spin-down evolution is sensitive to the specific propeller 
mechanism rather the initial spin period of the NS. 



2. Spin-down models of NSs in HMXBs 

2.1. The ejector phase 

Normally a newb orn NS in a bi nary first appears as a radio pulsar (or ejector) after the 
supernova explosion ( ]Lipunovlll992l ). In this phase the spin-down of the NS is due to the 
loss of its rotational energy dominated by magneto-dipole radiation and the outgoing flux 
of relativistic particles. If the NS's companion is a high-mass star, the stellar wind matter 
from the companion within the gravitational radius of Rq = 2GM/V 2 will be captured by 
the NS at a rate M = nR^pV . Here G is the gravitational constant, M the mass of the NS, 
p the density of the wind at Rq, and V the velocity of the NS relative to the stellar wind, i.e. 
^ = (K>rb+Kv) 1 i where V or ^ and V w are the orbital velocity of the NS and the wind velocity, 
respectively. The pressure of the outgoing radiation and particles is larger than that of the 
i ncoming matter a t Rq. The energy loss rate can be expressed as E = — /i 2 f2 4 (l + sin 2 a)/c 3 
( iSpitkovskyl 120061 ) . where p = BR 3 is the magnetic moment of the NS (B and R are the 
surface magnetic field and radius of the NS, respectively), fl the angular velocity, and a the 
inclination angle between the magnetic and rotational axes. Thus the spin-down rate in the 
ejector phase is 

■ _ An 2 B 2 R 6 (l + sm 2 a) 

IP? ' () 

where / is moment of inertia of the NS. As the NS spins down and the outgoing pressure goes 
down, the transition to the supersonic propeller phase will occur when the two pressures are 
in balance. The spin period of the NS at the transition point is 
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B 2 R 6 (1 + sin 2 a) 
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(2) 
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2.2. The propeller phase 



Once the wind matter crosses the gravitational radius Rq the propeller phase starts. If 
the plasma enters the light cylinderical radius the pulsar mechanism will switch off and the 
incoming matter will form a quasi-static atmosphere surrounding the NS. At this moment 
accretion does not occur since the magnetosphere radius R m = (/z 2 / ' M\j2GM) 2 / 7 is larger 
than the corotation radius R co = (GM/Q 2 ) 1 / 3 of the NS. The infalling material is stopped at 
the magnetosphere by the centrifugal barrier, which prevents material from accreting onto 
the NS. The ejected material will carry away the angular mome ntum of the NS and de c elerat e 
its spin. This so-called propeller effect was first introduced by llllarionov &; Sunyaevl ( 119751 ). 



Davies et al.l ( 119791 ) pointed out that, according to the Mach number Ai = QR m /c s (here 
c s ~ (GM/i? m ) 1//2 is the sound velocity at R m ) at the magnetosphere, the propeller phase 
can be subdivided into two cases: supersonic propeller and subsonic propeller. Accordingly 
the above mentioned propeller mechanism is related to the supersonic propeller since the 
Mach number Ai > 1. This phase ends when R co = R m (i.e., Ai = 1) and the corresponding 
spin period 

P eq = 2 11 / 1 V 6/7 M~ 3/7 (G'M)- 5 / 7 (3) 



is called the equilibrium period. Further works ( lArons fc Leal Il976l ; Eisner &: Lamb! 119771 ) 
showed that, unless the material outside the magnetosphere is able to cool, accretion is 
unlikely to happen. Thus, even with R co > R m the propeller stage will succeed once the 
energy deposition rate is larger than the energy loss rate of the surrounding shell, and keep 
removing angular momentum from the NS. Because the Mach number Ai < 1, this stage is 
called the subsonic propeller. This process will cease if the loss rate of the rotational energy 
can no longer support the surrounding atmosphere against cooling, then the atmosphere will 
collapse and the NS enters the accr etor stage. The spin period of the NS at this point is 
so-called the break period, given by (jDavies &: Pringldll98ll ; Ilkhsanovll200ll ) 



-Pbr 

where /x 30 = /u/10 30 Gcm 3 , and Mi 6 
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It should be noted that the supersonic propeller can occur in both wind-fed and disc-fed 
cases. However, there is no consensu s on the angular momentum loss rate of a NS during the 
propeller phase (Ipavies et al.lll979f). Here we adopt a general formulation of the spin-down 
torque as follows ( IMori fc Rudermanl 120031 ) . 



iCl = -MR 2 m tt K (R n 



(5) 



where 7 is a parameter ranging from —1 to 2, and its value reflects various propeller mech- 
anisms and spin-down efficiencies. For the supersonic propeller, 7 =—1, and 1. When 
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7 = -1, the matter is assumed to be ejec ted with the escape velocity at R m , i.e., v esc (R 

ml 

y^2GM/R m (llllarionov fc Sunyaevlll975l ). and the spin-down torque is calculated based on 
energy budget. The energy loss rate is IQQ = — (l/2)Mt> 2 sc (P m ) = — M[R m Q K (R m )} 2 , where 
^k(-Rhi) is the Keplerian angular velocity at R m . When 7 = and 1, the matter is assumed 
to be ejected at the escape v elocity v esc (R m ) ( jDavidson fc Ostrikerlll973l ) and the rotating 
velocity R m Q (jShakuralll975l ) of the magnetosphere at P m , respectively, and the torque is 
derived under the angular momentum budget. The corresponding angular momentum loss 
rate is iCl = —MR m (2GM/R m ) 1 / 2 = -2^ 2 MR 2 m n K (R m ) and Itt = -MR^Vt, respectively. 
The value of 7 for the subsonic propeller phase is 2, in which the rotational energy of the NS 
is assumed to be dissipated at a rate of IVLVt = — Mi?^f2^(i? m )[n/f2 K (-R m )] 3 , and the result- 
ing torque is IQ = — MP^fiK^PmH^/^K^Pm)] 2 - Thus the spin-down rate in the propeller 
stage can be summarized as 



. (2ny-\GM)^MR m 

r 7P7-2 ' W 

As an illustration, we consider a 1.4 M NS with an initial spin period of 0.01 s, a surface 
magnetic field of 10 12 G, and an accreting rate of 10 16 gs _1 . The corresponding characteristic 
spin-periods are P e j ~ 0.24 s, P cq ~ 6.7 s and Pbr — 81.5 s respectively. The timescales in the 
supersonic propeller phase varies from 30kyr to 20Myr as 7 decreases from 1 to — 1, and 
in the subsonic propeller phase which has a spin-down rate irrelevant with P the spin-down 
ti mescale is ~ 90kyr. In th e above calculation we use 300kms _1 as the relative velocity (see 



1+37 



Raguzova fc Lipunovlll998l ). 



2.3. The accretor phase 



Steady wind accretion onto the NS starts at P > Pbr- In this phase the spin period could 
be further ch anged since the wind matter possesses some angular mo mentum. However both 



Matsuda et al 



1992 



observations (Bildsten et al.lll997h and numerical calcula tions (e.g.. iFryxell &: Taamlll988 



Azner fc Bornerlll995l ; lRuffertlll999l ) have shown that the efficiency of 



angular momentum transfer in wind accretion is quite low, with alternating short-term spin- 
up and spin-down. Thus one may expect that the present spin periods of wind-fed X-ray 
pulsars are not significantly different from the Pb r achieved earlier. 



Recently IShakura et al.l (120121 ) proposed a model of subsonic quasi-spherical accretion 
onto a slowly rotating NS in HMXBs with low X-ray luminosities (Lx < 10 36 ergs _1 ). In 
this model the accreting matter settles down subsonically onto the rotating magnetosphere, 
forming an extended quasi-static shell around it. The angular momentum can be removed 
from or injected into the NS depending on the sign of the specific angular momentum of the 
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falling matter. In the case of moderate coupling between the plasma and the magnetosphere, 
from the torques acted on the NS due to both the magnetosphere-plasma interaction and 
accretion, the changing rate of the spin period is derived to be (see also PT12) 



P 2 



2ixl 



Mj 6 3+2n)/n - GMT 



(7) 



where A ~ 2.2 x lO^P^Pe) 1 / 11 ^^ 300 , and C ~ 5.4 x lO 31 !^^^) 13 / 11 /^. 
Here B 12 = P/10 12 G, R 6 = P/10 6 cm, P 1000 = P/lOOOs, P orbi30 o = P OT b/300hr, V 300 = 
V/300 kms -1 . The constants K\ and n are set to be 40 and 2 respectively. It is seen that 
there is an critical accretion rate at which P = 0. 



3. Estimate of the magnetic field 

The NS magnetic field is a critical parameter in the spin-down models. Before investi- 
gating the spin history of SXP 1062, we need to know the information about its magnetic 
field strength. The cyclotron features in the X-ray spectra provide the most direct and ac- 
curate way to measure the magnetic field strengths of accreting NSs. Unfortunately, they 
have not been detected in SXP 1062. Nevertheless, there are several other ways to estimate 
the NS magnetic field from its spin period and period derivative, though model dependent. 

One of the hints comes from the young age of the SNR associated with the NS. This 
requires that the lifetime of the ejector phase (usually much longer than that of the propeller 
phase) must end within a few 10 4 years. Assuming that the magnetic field has changed little 
during this phase and that the initial spin period is much smaller than P e j, one can estimate 
the timescale of the ejector phase to be (PT12) 

Tcj = ~ 15 ^Vaoo 172 ^ 1 Myr. (8) 

This value is about two orders of magnitude larger than the estimated age of SXP 1062, 
unless B12 > 100. This means that SXP 1062 must have possessed very strong magnetic 
field. 

PT12 further assume d that the NS in SXP 1062 is spinning at the equilibrium period as 



described in the model of IShakura et al.l (120121 ). and derived the current magnetic field to be 
B±2 < 10 13 G using Eq. (7). Accordingly they suggest that the NS magnetic field must have 
been str onger in the pas t and then decayed to its present, normal value. It is noted that the 



model of IShakura et al.l (120121 ) has quite a few parameters whose magnitudes are uncertain. 



For example, the value of K%, which relates the poloidal (B p ) and toroidal (B^) magnetic 
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field components, is found to be ~ 40 in IShakura et al.l (120121 ). This will result in 3> B p 
during the accretor phase, an d it is not known wheth er the magnetic field configuration can 
remain stable in this case (cf. lAlyl Il985l : IWang Ill995l ). 



The extraordinary large spin-down rate of SXP 1062 can b e used to put useful constraint 



on the magne t ic field of the NS. As sh own by many authors (ILynden-Bell fc Pringldll974 



Lipunovlll982l ; lBisnovatyi-Kogan!ll99ll ). the maximum spin-down torque exerted on a NS in 
either disc or spherical accretion is 



m = -k 



R 3 '' 



(9) 



where k < 1. To account for the spin-down rate measured in SXP 1062, the NS magnetic 
field has to be 

B~3x 10 14 /«- 1/2 M 1 1/ 4 2 / 4 1 5 /2 J R 6 - 3 (P/100 syr- 1 ) 1 / 2 G, (10) 

where M 14 = M/1.4M , and J45 = I/10 45 gem 2 . The same result can be obtained if the 
spin-down torque in the subsonic propeller phase (jDavies fc Pringldll98ll ) is used. 



Another efficient spin-down mechanism was proposed by llllarionov fc Kompaneetd ( 119901 ) 
They argued that there could be outflows from the NS magnetosphere caused by heating 
of hard X-ray emission of the NS if the X-ray luminosity falls in the range of ~ 2 x 10 34 



ergs — 3 x 10 ergs . Compton scattering heats the accreted matter anisotropically, and 
some of the heated matter with a low density can flow up and form outflows to take the 
angular momentum away. The corresponding spin-down torque is 



Z7l 



ID 



Here M out is the mass outflow rate (no larger than the mass transfer rate) and x is the solid 
angle of the outflow. This gives the magnetic field to be 



B * 3.6 x 10 14 (g )- 7 /«/J/ 8 M{'^R«\^_ 



_ r 3/8 ( 



100 syr- 



1062 s' 



- ? / 8 G. (12) 



The above estimates show that SXP 1062 could be an accreting magnet ar. Similar 



conclu sions have also been drawn for other X-ray pulsars in HMXBs. iDoroshenko et al. 
( I2010l ) reported the spin history of the 685 s X-ray pulsar GX 301—2, and found it spinning 
down at a rate v ~ 10~ 13 Hzs" 1 . Reig et al. (2012) showed that the measurements of the 
spin period (5560 s) of 4U2206+54 imply a spin-down rate of v ~ — 1.5(±0.2) x 10 -14 Hzs -1 . 
Using the above spin-down mechanisms to explain the spin-down rat es also leads to very 
strong magnetic fields (> 10 14 G) in these NSs (see also lLipunovlll982l ). 
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Ikhsanov fe Fingerl (120121 ) suggested an alternative interpretation for the rapid spin- 
down in GX 301—2. They showed that if the accreting material is magnetized, the magnetic 
pressure in the accretion flow increases more rapidly than its ram pressure, and under certain 
conditions the magnetospheric radius 



1.5 x itfofJML 1 3 Rf /13 T 6 - 2/13 Mtt 3 M 16 - -cm, 



-4/13 



(13) 



is considerably smaller than the traditional magnetospheric radius. Here a = O.lao.i is the 
efficiency parameter of Bohm diffusion, and T = 10 6 T6 K is the plasma temperature at the 
magnetospheric boundary. The spin-down torque applied to the NS is found to be 



K m fM 



( -^co Rmca, ) 3 ^ 



(14) 



where K m is a dimensionless efficiency parameter for the magnetic viscosity coefficient, and 
< K m < 1. The above equation can explain the spin-down of GX 301—2 with a normal 
field of a few 10 12 G if n m ~ 0.1. In the case of SXP 1062, it results in the estimate of the 
magnetic field to be 



i?^2xio i v :r /i 7 4 i f i7 M 1 



13/17 r 13/17 ,,8/17 r 13/17 D -3 s >-6/17^,-3/17. 



P 



l 45 



l 16 



-6 



100 syr- 1 



,13/17/ 



P 



1062 s' 



- 13/17 G, (15) 



where «o.i = w m /0.1. In the same way, llkhsanovl ( 120121 ) estimated the magnetic field of SXP 
1062 to be ~ 4 x 10 13 G by assuming K m — 1. Even this limiting value is comparable to the 
quantum critical field Bq = 4.4 x 10 13 G. 

According to the above arguments, in the following we assume that the current magnetic 
field of SXP 1062 is > 10 14 G. As to the evolution of the magnetic field, we consider two 
kinds of models. First we assume that the magnetic fi eld was initially stron ger and adopt a 
phenomenological model for the magnetic field decay (iDall'Osso et al.l 120121 ) 



dB 



-AB 



l+ct 



B 



(16) 



where the field decay timescale t&{B) = (AB a ) 1 , and A and a are constants. The solution 
of the above equation in the case of a ^ is 



B = Bi(l + at/r^y 



-l/a 



(17) 



where B\ is the initial field strength and 7a,; = (AB?)~ l . iDalFOsso et al.l (120121 ) show that, 
to be compatible with the observations of magnetar candidates, the magnetic field should 
decay on a timescale of ~ 10 3 yr for B ~ 10 15 G, with a decay index most likely within 
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the range 1.5 < a < 1.8. Here we adopt the initial magnetic field as 7 x 10 14 , 3 x 10 14 and 
10 14 G, with a = 1.6 and r d)i = 10 3 /B? 15 yr where B i>15 = B-J10 15 G. 



On the other hand, the observe d braking indices for several young radio pulsars have 



been measured and are all less than 3 ( 



20051 . 120061 : iLivingstone fc Kaspill2011 



vne et alil993 



1996 



Weltevrede et al. 



Kaspi et al.lll994j ; ILivingstone et al 



201l[ ) , suggesting that the NS mag- 
netic fields may be increasing. In particular, the brakin g index of the high-fie ld (5 x 10 13 G) 



radio pulsar PSR J1734— 3333 was found to be 0.9±0.2 (jEspinoza et al.ll201ll ). implying that 
this pulsar may soon have the rotational properties of a magnetar. In the second approach, 
we adopt a field growth model in the following form 



B = 

with Bi = 3x 10 12 G, r = 10 3 yr, and a 
at t = 10 4 and 4 x 10 4 yr, respectively. 



5i(l + Vr)°, 
= 1.45, so that B 



8.5 x 10 13 G and 6.3 x 10 14 G 



In Figure 1 we show the model evolution of the magnetic fields. 



4. Spin evolution 



A newborn NS is usually rotating rapidly. However, lHaberl et al.l (120121 ) suggested that 
SXP 1062 could have been born with a period much larger than 0.01 s. Some central compact 
objects (CCOs) in supernova remnants which have spin periods ra nging from ~ 0.1 to ~ 0.5 s 
( Zavlin et aliboOO ; Gotthelf et al. 2005 ; Gotthelf fc Halpern 2009 ) seem to support this point 
of view, since there is evidence that the spin periods of these sources are very close to the 
initial ones. Thus in our model we take 0.01 s, 0.5 s and 6.5 s as the initial period of the NS in 
order to examine whether it can significantly influence the spin-down evolution. We use the 
ultra-long initial spin period of 6.5 s because this value is larger than P e j with B = 7 x 10 14 G 
in the ejector phase, so that the NS will directly enter the supersonic propeller phase after 
the SN event. 



It was shown by lArons &: Leal (119761 ) and lElsner fc Lambl (119771 ) that, for stable accre- 
tion to occur, the plasma at the base of the NS magnetosphere should become sufficiently 
cool, so that the magnetospheric boundary is unstable with respect to interchange instabil- 
ities. This can be realized only if the spin period of the star exceeds the break period P^, 
and the X-ray luminosity is larger than 



3 x 10 



^X^e^ergs- (19) 

If B > 10 14 G, SXP 1062 should be in the subsonic propeller phase. However, it is not 
clear whether the picture of the subsonic propeller can be applied to BeXBs, due to the 
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following reasons. (1) The mass accretion in BeXBs is now believed to be triggered by 
Roche- lobe overflow of the Be discs which is trun c ated by th e NS through a tidal torque 
( lOkazaki fe Negueruelal l200ll ; lOkazaki et al.l 120021 ; iReigl 1201 ll ). thus is deviated from the 
traditional Bondi accretion in supergiant HMXBs. This means that the NS in SXP 1062 is 
probably surrounded by a (quasi-)disc rather a quasi-static, spherical atmosphere. (2) Even 
in the spherical wind-fed case, the spin period - orbital period correlation in BeXBs seems 
to be well accounted for by assuming that the NSs are spinning at the equilibrium periods 
described by Eq. (3) ( Corbetlll984j ; IWaters &: van Kerkwiilj[l98M . Thus in our calculations 
we don't consider the subsonic propeller phase, and assume that the evolutionary sequence of 
the NS is ejector- supersonic propeller-accretor. We use different 7 to calculate the spin-down 
torque in the supersonic propeller phase, and take the equilibrium period (Eq. [3]) as the 
final period (i.e., the period does not change during the accretor phase unless the magnetic 
field changes). 

In Figures 2 — 4 we show the calculated results corresponding to different initial spin 
periods of the NS. Here we take the NS mass to be M = 1.4M , the inclination angle a = 90°, 
/ = 10 45 g s -1 , and R = 10 6 cm. The relative wind velocity V is set to be 300 km s -1 , and the 
accretion rate is fixed to be 10 16 gs _1 . The three thin lines (from top to bottom) describe the 
spin evolution with initial magnetic field of 7 x 10 14 G, 3 x 10 14 G, and 10 14 G undergoing field 
decay, respectively; the thick line is for the field growth model with initial field of 3 x 10 12 
G. The solid, dashed, and dotted lines represent the ejector, propeller, and accretor phases, 
respectively. 

We notice that the time spent in the supersonic propeller phase is sensitive to the value 
of 7 which reflect different spin-down mechanisms in the supersonic propeller phase as we 
mentioned before. In the case of 7 = — 1, where the spin-down torque is most inefficient, the 
NS has not evolved out of the supersonic propeller phase at the age of the SNR, even with 
a superstrong magnetic field. Our results are not sensitive to the initial spin period of the 
NS. Thus even for the case of ultra-long initial spin period there is no significant change in 
the final NS period. 



In other cases the NS can successfully reach P eq when t = 10 — 40 kyr. Since P ( 



cq 



depends on the magnetic field, we can see that the magnetic field determines the final spin 
period which the NS can achieve, and the value of 7 determines the evolutionary timescale. 
In the field decay model, the spin period remains invariant once it reaches P eq , since we 



1 Additionally, population synthesis calculations by iDai fe Lil (|2006l ) showed that the distribution of the 
spin and orbital periods of X-ray pulsars in supergiant HMXBs can be ro ughly explained w ithout requiring 
that an X-ray pulsar emerges after the subsonic propeller phase (see also lStella et al.lll986l ). 
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assume that during the accretor phase the long-term, net torque from the wind is small. In 
the field growth model, the spin period keeps increasing with P eq in the final stage, since 
the increase of B always breaks the instaneous equilibrium and causes a spin-down torque. 
It is seen that B > 10 14 G can fulfill the requirement of SXP 1062 in both models. This 
result favors that in the supersonic propeller phase matter is ejected at R m with the escape 
vel ocity under angular mome ntum conservation, consistent with the numerical calculation 
by IWang fc Robertson! (119851 ) . 



5. Discussion and conclusion 



The newly found Be/X-ray binary SXP 1062 is believed to be the first X-ray pulsar 
associated with a SNR, which shows a combination of young age and long spin period that 
cann ot be explained by a typical NS. Previous studies (IHaberl et al.ll2012t iPopov fc Turolla 
20121 ) explored its possible origin invoking initially long spin period or ultra-strong magnetic 
field. Here we discuss the possibility that SXP 1062 is an accreting magnetar with B > 10 14 
G, and examine in this case how the properties of the NS (i.e. initial spin period, magnetic 
field and its evolution) and the spin-down torques can be constrained. 



Other candidates of accreting mag net ars in binaries includ e 4U 2206+54 ( IFinger et al. 



20101 ; iReig et al.ll2012l ) and GX 301-2 f lDoroshenko et a 



201ph . However, it is controversial 



whether they really possess ultra-strong magnetic fields. IWangl (120101 ) reported the existence 
of two cyclotron absorption lines at ~ 30 and 60 keV in 4U 2206+54, and derived a magnetic 
field of 3.3 x 10 1 2 G, although no s i gn of this feature has been detected in other observations. 
Observations of lLa Barbera et al.l (120051 ) showed the cyclotron resonanc e scattering feature 
at ~ 35 — 45 keV in GX 301—2, suggesting the field strength of 4x 10 12 G. llkhsanov fc Finger 



(120121 ) proposed a magnetic wind accretion model for GX 301— 2 to account for the difference 
in the field strengths derived from the cyclotron lines and from the spin-down rates. In the 
case of SXP 10 62, it is f o und t hat the magnetic field may be at least strong as ~ Bn 



in the model of llkhsanovl (120121 ). For a dipole magnetic field of ~ 10 14 G, the electron 
cyclotron line would appear at E > 1 MeV, but a proton cyclotron line would appear at E ~ 
0.5(5/10 14 G) = 0.3 keV. Although a line with this energy should be obser vable with XMM - 
Newton detectors, it is in a region affected by strong interstellar absorption (IReig et al.ll2012l ). 
Currently, no sign ificant lines have been detected in the persistent emission of magnetars 
(lMereghettill2008h . 



If SXP 1062 is or has been a magnetar, the association between the SNR and SXP 
1062 may provide an opportunity to investigate the formation and evolution of magnetars. 
Vink fc Kuiperl (120061 ) showed that there is no evidence that magnetars are formed from 
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rapidly rotating proto-neutron st ars. The SNR associated with SXP 1062 is one of the 



faintest SNRs known in the SMC (IFilipovic et al.ll2005l 120081 : lOwen et aUEOllh . This seems 
to in line with the finding of IVink fc Kuiper that their formation may not be accom- 

panied with extraordinarily bright supernovae. However, it is known that the brightness of 
SNRs depends strongly on the density of the environment. Nevertheless, the age of the SNR 
can be used to set useful constraints on the timescale of magnetic field evolution, either due 
to field decay or growth. 

Since the long spin period is most likely to be reached during the propeller phase, 
the age of the SNR also plays a role in testing the efficiency of the spin-down torques in 
different propeller mechanisms. Our results seem to rule out the model with 7 = —1, and 
prefer larger values of 7 which correspond to more efficient prope ller spin-down. Recent 
2- a nd 3-dimensional magne tohydrodynamic (MHD) simulations by iRomanova et al.l (120051 ) 
and Ustvugova et ali ( 2006 ) on disc-fed NSs suggest tl oc — fi 2 for propeller-driven outflows. 
Toropina et al.l (120101 ) investigate the spinning-down of magnetars rotating in the propeller 
regime with axisymmetric MHD simulations, and find Q oc — Q 1 ' 5 . It should be noted that 
the mass transfer rate is assumed to be constant throughout our calculations, but in reality 
it must have varied with the orbital motion of the NS. For instance, an eccentric orbit 
may result in alternation among the ejector, propeller and accretor phases. The sporadic 
outburst behavior will further complicate the spin-down evolution of the NS. This means 
that the calculated evolutionary sequence in our model and the values of 7 should be taken 
as an illustration and lower limits, respectively. However, both the high spin-down rate 
and the young age of SXP 1062 provide strong evidence that the binary indeed harbors or 
harbored a magnetar, and a effective spin-down mechanism is required. We expect further 
observations to confirm the long-term spin behavior of SXP 1062. 
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Fig. 1. — Evolution of the NS magnetic fields in different models. 
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Fig. 2. — The spin-down evolution for a NS with a 0.01 s initial spin period. In this and 
following figures, solid, dashed and dotted lines represent the ejector, supersonic propeller 
and accretor phases, respectively. Symbol * represents the transition of evolutionary stages. 
The shaded area and the horizontal line mark the age and the spin period of SXP 1062 
respectively. From the top to the bottom, the initial magnetic field is B\ = 7 x 10 14 G, 
3 x 10 14 G, 10 14 G, and 3 x 10 12 G. 
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Fig. 3. — The spin-down evolution for a NS with a 0.5 s initial spin period. Other parameters 
are similar to in Figure 2. 
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Fig. 4. — The spin-down evolution for a NS with a 6.5 s initial spin period. Other parameters 
are similar to in Figure 2. 



